Contour dynamics model for electric discharges 
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We present an effective contour model for electrical discharges deduced as the asymptotic limit of 
the minimal streamer model for the propagation of electric discharges, in the limit of small electron 
diffusion. The incorporation of curvature effects to the velocity propagation and not to the boundary 
conditions is a new feature and makes it different from the classical Laplacian growth models. For 
the first time, the dispersion relation for a non planar 2-D discharge is calculated. The development 
and propagation of finger-like patterns are studied and their main features quantified. 

PACS numbers: 51.50.+V, 52.80.-s 
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The electrical breakdown of various media is usually 
caused by the appearance and propagation of ionization 
waves and, in particular, streamers. The difficulties for 
the understanding of this important and interesting phe- 
nomena are both experimental and theoretical. From 
experimental point of view, the historical lack of quan- 
titative data, from small to large scales discharges, are 
mainly due to the fast time scales of the processes in- 
volved. Recently this situation is starting to change with 
the development of new experimental techniques. On the 
other hand, the theory is still far from being complete. 
Even in the simplest hydrodynamic approximation for 
describing the phenomena, it is necessary to solve a non- 
linear system of balance equations, together with Pois- 
son equation, which poses a challenging problem both 
numerical and analytical. Among the recent progress in 
the understanding of the propagation mechanism we can 
mention: the study of stationary plane ionization waves 
[U, Q > nrs t self-similar solutions for ionization waves in 
cylindrical and spherical geometries [3, [j| , and the for- 
mation of streamers as the result of the instability of 
planar ionization fronts @-@|- in the simplest hydrody- 
namic approximation, the fronts are subject to both sta- 
bilizing forces due to diffusion which tend to dampen out 
any disturbances, and destabilizing forces due to electric 
field which promote them. 

In this letter we introduce and justify theoretically a 
contour dynamical model which allows us to make some 
progress in 2-D and 3-D more general situations. Here 
we will use ideas coming from the context of electrohy- 
drodynamics. More precisely, we will arrive at a model 
similar to the celebrated Taylor-Melcher leaky dielectric 
model for electrolyte solutions [9j, but suitably adapted 
to the context of electric (plasma) discharges. We will ob- 
tain the dispersion relation for 2-D discharges like those 
described in some experiments 

For our contour model, in the asymptotic limit of small 
diffusion D -C 1 (electron diffusions are typically of order 
0.1 m 2 /s and, since typical velocities and streamer size 
[1] are 10 5 m/2 and 10~ 2 m we get D of order 10 -3 ), the 



interface will move according to 
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being the normal velocity and k twice the mean cur- 
vature and a(E+) to be defined below. The free charge 
density on the boundary will be concentrated in a dif- 
fuse layer, and for the limit considered, its behavior will 
be characterized by a negative charge surface density a 
given by 
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where j~ is the current density coming from the ionized 
region (Q) to its boundary (dfl) in the normal direction 
and g is related to the resistivity of the plasma. E„ is 
the normal component of the electric field at the inter- 
face when approaching it from the region without /with 
plasma, i.e 
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with V found by solving 

AV = aS(r-r dn ). (4) 

Notice that in the case g^ 1 1, we arrive to Lozansky- 
Firsov model [l(| with a correction due to electron dif- 
fusion, meanwhile in the limit D — we arrive at the 
classical Hele-Shaw model. Such a model is known to 
posses solutions that develop singularities in the form of 
cusps in finite time [ll| but, when regularized by sur- 
face tension corrections, the interface may develop yari- 
ous patterns including some of fractal-type (see [12| for 
a recent development and references therein). 

Let's derive the equations (pj and ([2]). We start from 
a minimal model of the discharge that in dimensionless 
units reads 

1 ■ (n e E + D Vn e ) 
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n e a(|E|), 
n e a(|E|), 
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witha(|E|) = |E| exp(— 1/|E|). Next we derive a contour 
dynamics equation for the evolution of the interface be- 
tween the region occupied by the plasma and the region 
free (or with a very small density) of plasma. Note that 
this interface does not exist as such but, as in the case 
of planar fronts, the region where n e jumps from values 
close to 1 to values close to has a thickness O(D^) 
Here we will give an equation, valid in the asymptotic 
limit D < 1. 

We take a level surface of n e representing the interface, 
and introduce local coordinates r (along the level surface 
of n e ) and v (orthogonal to the level surface of n e ). Since 
the front thickness is O(D^) it is natural to introduce 
v = D^Xi so we can write approximately 

D ^^^4 x +D ( A --^m) +0iDi) 

where A_l is the transverse laplacian and k is twice the 
mean curvature in 3-D or just the curvature in 2-D (see 
e.g. [HI]). Hence from ([5]) and (0 
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We expand now |E|and n p at the diffusive boundary 
layer as explained in Q and Q, so it is straightforward 
to end with 

dn e _ dn e ( E v „ , — — ; i— \ dn e _ . „ i . 



.(8) 

where a curvature term appears and is relevant provided 
1 <C k <C D~*. Eq. ([8]) is a transport equation for the 
electron density so that the level line of n e that we have 
chosen to describe the evolution of the interface moves 
with a normal velocity given by (UJ . Notice that the tan- 
gential velocity does not change the geometry of the curve 
during its evolution. Nevertheless, tangential exchanges 
of charge affect the evolution through the dependence of 
vn on E v . 

Next we describe the charge transport along the inter- 
face. We trace a "pillbox" T> around a surface element. 
The region T> will be such that At ^> Ai/ and will con- 
tain both the diffusive layer for n e and the region where 
n e — rip (the net negative charge) has significant values. 
If we subtract © from ([5]) and integrate over V we find 
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where we have used that n e — > for x — v I^TD 3> 
1, |Vn e | — > for |x| ~> 1 and we have neglected the 
contributions of the lateral transport of charge by E T in 
comparison with exchange of charge with the bulk by E v . 



This assumption is also common in the Taylor-Melcher 
model and we will follow it. Notice that 
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and the length (or area) element of the surface will suffer 
a change during interface deformation given by 
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We can also include a source (an insulated wire inside 
the plasma, for instance at Xo , carrying a electric current 
I(t)). This source will create a current density inside the 
plasma and as quasineutrality is fulfilled, we will have for 
the interior of fi 



V-j = /(i)<5(x-x ). 



(12) 



By putting ©-(HU) together we can finally write the 
equation (2}, where g is the effective resistivity of the 
medium, given in this case as g^ 1 = lim x= _ oc n e and 
E~ is the normal component of the electric field when 
approaching the interface from inside the plasma region. 
j~ is normal component of the current density given by 
(TT2")) arriving at the interface. There is a jump in the nor- 
mal component of the electric field across the interface 
given by 



E, 



E~ = -a. 
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Eq. @ will provide the surface charge density a as a 
function of time. From it, we can compute the electric 
field and move the interface with (fT]). Two limits can be 
easily identified in the case of j « 0: a) the limit of large 
conductivity 

g- 1 > 1, E~ =0^> V = C(t). 

so that the interface is equipotential and b) the limit of 
small conductivity 



g- 1 « 1, ^ {a At) = 



a At = Const. 



where the charge contained by a surface element is con- 
stant and the density only changes through deformation 
(with change of area) of the interface. 

Next we study as an application a 2-D case. The 2-D 
case is not at all academic as there have been experiments 
where a 2 dimensional streamer discharge is created on 
a surface and branched pattern are observed. In Fig. [1] 
we can see a numerical simulation intended to mimic the 
experiment of [l[ using dTJ and ([2} . This set of equations 
also allow us to calculate the dispersion relation for the 
front stability and give some analytical insight. 

Let's write the position and the charge "surface" den- 
sity of the interface as 



r{6,t)=r(t)+6S(9,t) 
Q{t) 
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FIG. 1: The evolution of an initial arbitrary front with a 
constant net charge Q = 30 and D = 0.1. 



where r(t) is the solution of the equations for the radial 
symmetrical front, Q(t) — J I{t)dt and 8 is our book- 
keeping small parameter for the expansions that follow. 

We will start calculating the correction to the electric 
field due to the geometrical perturbation of the surface 
and the extra charge deposited on it. The electric poten- 
tial will be up to delta order $(x) = V(x) + 8V p (x) being 
V the solution for the symmetrical problem. The term 
V p (x.) will satisfy the equation AV P = 0(8) and further 
on, we will change coordinates to x — > x = xr(t)/r(0, t) 
so the perturbed surface becomes a disk of radius r(t). 
Now at zero order V p satisfies the Laplace equation with 
the boundary being a disk. Hence we have in the new 
polar coordinates 



V P (f, 




(16) 



where we have imposed the conditions of V p remains finite 
at the origin and becomes zero at very large distances. 
The potential at the surface position coming from the 



exterior is $(x^ 



-Q(t) log(r(0,t))/27r + <%O 



and 

from the interior is 3>(x~) = C(r(t),t) + i5V^(xj) where 
C(r(t),t) is a function independent of 6. So if we write 
the surface perturbation as 



S = s n (t) cos(n(9), 



(17) 



n=l 



imposing the condition of continuity for the potential, 
we find that the coefficients of the potential in (fT6"| are 
related by 



, i Q(t) 

Ipn = + T, s r. 
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Now we can calculate the electric field to 8 order. Chang- 
ing back coordinates, from x to x, the normal compo- 



nents of the electric field at both sides of the surface are 
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Then, the jump condition (fT3|) together with (|T5|) give 
the following expression for the charge perturbation 
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Q(t) 

2nr 



cos(n6»), (19) 



where r = r(t) is the evolving radius of the unperturbed 
circle, i.e. the radially symmetric solution to (H},©. 

We still need the expression of the curvature at 8 or- 
der to find the dynamics of the front. In polar coordi- 
nates, after introducing the perturbation (|14[) . the cur- 
vature turns out to be 



-8 + 0(S 2 



_ r 2 +2rS8-rS ee 8 + 0(8 2 ) _ 1 Sj_ 
(r 2 + 2rS 9 S + 0(8 2 ))? r ^ 
and the normal component of the velocity 

dr(t) F dS(0,t) 

Now Eqs. (HJ) and @ after introducing the perturbations 
to first order, defining e = D, and doing some algebraic 
manipulations yield 
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Thus the time evolution of each particular mode has been 
obtained and it is governed by (|20|) and (f21~j). First let's 
study the limit of ideal conductivity. It makes g — > 0, so 
we can see from ([2~Tj) that cp n 0. Physically it means 
that in the limit of very high conductivity, the electric 
field inside goes to zero (E~ — > 0), as we approach to 
the behaviour of a perfect conductor. If we consider that 
Q(t) = Qo is constant or its variation in time is small 
compared with the evolution of the modes (which also 
implies I it) — > 0), and the same for the radius of the front 
r(t) = Tq, we can try a solution s n = exp(w n t), ip n = 0, 
to our system, and we get a discrete dispersion relation 
of the form 
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Next we consider the limit of finite resistivity, but such 
that the total charge is constant at the surface, or varies 
very slowly. Writing (f2"Tj) as 
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FIG. 2: Dispersion relation for the discrete modes of a per- 
turbation keeping the charge constant for different charges 
Q — —Qo- In the inset, the case Q = 10 is compared with the 
case of ideal conductivity (cross points). 



we have now 



dt 



Qo ds n _ J_ 
Atttq dt 2g 



For a small enough conductivity, so no extra charge 
reaches the surface, we find ip n ■ 
s n = cxp(w„i), (|20| yields 



Qo 



and with 
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In a curved geometry we can see that the modes are 
discrete. If we compare (|22l) and (|2"3")l . for n = 1 cor- 
responds the same uj\ = which reveals the translation 
invariance of the front. However, for small n, there is a 
1/2 factor discrepancy in the dispersion curve between 
the two limits. In [14j . the origin of this pref actor was 
discussed for planar fronts (the dispersion relation for 
planar fronts was obtained in the case of constant charge 



in Q). So imposing constant charge at the surface, this 
factor is independent of the planar or curve geometry. On 
the other hand, imposing constant potential at the sur- 
face gives a factor of 1. The intermediate situations can 
be studied by solving the system ([20]) and (12T1) . Another 
consequence is that in both cases the maximum growth 
correspond to a perturbation with 



n cx \Q \/D, 



(24) 



implying that the number of fingers increases with the 
net charge and decreases with electron diffusion. 

To test the analytical predictions, we have calculated 
numerically the dispersion relation curves for these two 
cases. In Fig. [2] for different charges u)\ — 0, the slope 
increases with the amount of charge, the maxima moves 
to higher modes, and for larger n's the dispersion curves 
become negative as predicted by (|2"3")l . In the inset of 
Fig. [2] the dispersion relation curve is calculated for the 
case of constant charge and constant potential. The slope 
around the origin n = 1 is bigger for the case of ideal 
conductivity, when the interface is equipotential. The 
ratio of the slopes can be estimated as « 0.49 which is in 
agreement with our prediction of 1/2 (see [13]). 

To conclude, we have introduced a contour dynam- 
ics model a la Taylor-Melcher [9| for the streamers dis- 
charges. Our model contains as a particular limit the 
Lozansky-Firsov model [l(| with a correction due to elec- 
tron diffusion which effectively acts as a surface tension. 
In the limit D = the classical Hele-Shaw model [TTI.[l2l| 
is recovered. In this framework we have studied the sta- 
bility of 2-D discharges. It is our hope that this model can 
open the door for the study of 3-D problems of streamer 
discharges in realistic geometries. An extra advantage of 
this model is its relatively low computational cost and 
the high analytical insight which provides. 
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